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ABSTRACT 

Multi-epoch X-ray spectroscopy (0.3-25 keV) of the Seyfert 1 .2 galaxy Mrk 79 (UGC 3973) 
spanning nearly eight years and a factor of three in broadband flux are analysed. The data are 
obtained at seven epochs with either XMM-Newton or Suzaku. Comparison with contempo- 
raneous RXTE monitoring indicate that all flux states of Mrk 79 are represented by the data. 
The spectra are fitted in a self-consistent manner adopting a power law and ionised reflection 
to describe the broadband continuum. Modification of the spectra by a distant photoionised 
medium, seen predominantly in emission, are also included. Under the assumption that the 
inner disk is at the innermost stable circular orbit, our blurred reflection models give a spin of 
a = 0.7 ± 0.1. The reflection component in each spectrum is weaker than predicted by simple 
reflection models. If the illuminating X-ray emission is produced by flares above the disk that 
move at mildly relativistic velocities, however, diminished reflection is expected. Light bend- 
ing due to strong gravity near black holes can influence how the illuminating and reflected 
flux are observed; variations in Mrk 79 do not suggest that light bending is important in this 
source. 

Key words: galaxies: active - galaxies: nuclei - galaxies: individual: Mrk 79 (UGC 3973) - 
X-ray: galaxies 



1 INTRODUCTION 

The X-ray spectra of type I active galactic nuclei (AGNs) can ex- 
hibit significant variability from one epoch to the next and the gen- 
eral spectral shape can normally be associated with the flux state of 
the AGN. For example, in the low-flux state the power law emis- 
sion above ~ 2 keV is flatter and the excess above the power law 
at energies < 2keV (i.e. the so called soft-excess) is relatively 
stronger than in the high-flux state. In addition, the 2 — lOkeV 
band often displays complexity in the form of spectral curvature or 
spectral drops in the low-flux state (e.g. Gallo et al. 2004; Grupe 
et al. 2008; Bachev et al. 2009). However, defining the flux "state" 
for AGN is more complicated than for black hole binaries because 
the characteristic time scales are much longer in AGN. Long-term 
monitoring (i.e. years) is generally required to properly determine 
flux states. Gallo (2006) even suggested that the hypothetical UV- 
to-X-ray slope (ot ox ) may be better suited for defining the X-ray 
low state rather than X-ray flux alone. Consequently, it is challeng- 



ing to gain understanding of a particular system from single-epoch 
observation without information of its typical variability behaviour. 

An advantage of multi-epoch modelling is that model param- 
eters that vary on different characteristic time scales can be exam- 
ined and constrained accordingly. While black hole spin, disc in- 
clination and elemental abundances are not expected to change; 
power law spectral slopes can vary on day-long scales (perhaps 
more rapidly) (e.g. Sobolewska & Papadakis 2009; Iwasawa et al. 
2010); and the inner disc ionisation parameter and the absorber cov- 
ering fraction have been suggested to vary rapidly (10 3-4 s) (Ponti 
et al. 2006, 2010; Risaliti et al. 2009). 

In this work multi-epoch X-ray observations of the Seyfert 1.2 
galaxy Mrk 79 (UGC 3973; z = 0.022189) are examined to study 
the nature of its X-ray emission and variability. Mrk 79 is well- 
studied in the optical, benefiting from efforts to measure its black 
hole mass by reverberation mapping (e.g. Peterson et al. 2004). 
However, despite being relatively bright (> 10~ n erg cm -2 s _1 
between 0.5 — lOkeV) Mrk 79 is poorly studied in the X-rays. 
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Snap-shot observations of Mrk 79 were done in 2000 and 2001 with 
XMM-Newton, marking the first time the spectrum above 2keV 
was analysed (Gallo et al. 2005; 2006). The observations revealed 
a soft-excess, possible warm absorption, and a slightly broadened 
Fe Ka emission line. Gallo et al. (2005) report a possible 8 keV 
emission feature in the 2001 observation, but this is not confirmed 
with the deeper observations presented here. 

In the following section the observations and data reduction 
are described. In Section 3, the data sets are examined and charac- 
terised. The RGS spectrum is analysed in Section 4 to understand 
the nature of the ionised plasma. The high-energy spectrum above 
3 keV is examined in Section 5. A self-consistent model simultane- 
ously describing all the data is presented in Section 6. We discuss 
our results and conclusion in Section 7 and 8, respectively. 



2 OBSERVATIONS AND DATA REDUCTION 

Mrk 79 was observed with XMM-Newton (Jansen et al. 2001) on 
five occasions between 2000 and 2008, and once in 2007 with 
Suzaku (Mitsuda et al. 2007). A summary of the observations is 
provided in Table[TJalong with the nomenclature we adopt through- 
out this work. 

The EPIC pn (Struder et al. 2001) and MOS (MOS1 and 
MOS2; Turner et al. 2001) cameras were operated in small-window 
mode with the medium filter in place for all the observations except 
XMM6 when the thin filter was used. The Reflection Grating Spec- 
trometers (RGS1 and RGS 2; den Herder et al. 2001) also collected 
data during this time, as did the Optical Monitor (OM; Mason et al. 
2001). 

The XMM-Newton Observation Data Files (ODFs) from all 
observations were processed to produce calibrated event lists us- 
ing the XMM-Newton Science Analysis System (SAS v9 . . 0). 
Unwanted hot, dead, or flickering pixels were removed as were 
events due to electronic noise. Event energies were corrected for 
charge- transfer inefficiencies. EPIC response matrices were gen- 
erated using the SAS tasks ARFGEN and RMFGEN. Light curves 
were extracted from these event lists to search for periods of high 
background flaring. Significant flaring was detected during the fi- 
nal ~ 18 ks of XMM6 and data during these periods have been ne- 
glected. The total amount of good pn exposure is listed in TableQ] 
Source photons were extracted from a circular region 35" across 
and centred on the source. The background photons were extracted 
from an off-source region on the same CCD. Single and dou- 
ble events were selected for the pn detector, and single-quadruple 
events were selected for the MOS. The MOS and pn data at each 
epoch are compared for consistency and determined to be in agree- 
ment within known uncertainties (Guainazzi et al. 2010). For sim- 
plicity in analysis and presentation only the pn data are discussed 
here. 

The RGS spectra were extracted using the SAS task RGSPROC 
and response matrices were generated using RGSRMFGEN. The OM 
operated in imaging mode with various filters in place at different 
epochs. However, an ultraviolet image was obtained at every epoch 
in either the UV Wl (2450 - 3200 A) or UVW2 (1800 - 2250 A) 
filter. 

Suzaku observed Mrk 79 on 2007 April 03. The two front- 
illuminated (FI) CCDs (XIS0 and XIS3), the back-illuminated (BI) 
CCD (XIS1), and the HXD-PIN all functioned normally and col- 
lected data during this period. The target was observed in the XIS- 
nominal position. 

Cleaned event files from version 2 processed data were used 



in the analysis and data products were extracted using xselect. 
For each XIS chip, source counts were extracted from a 4' circular 
region centred on the target. Background counts were taken from 
surrounding regions on the chip. Response files (rmf and arf) were 
generated using xisrmf gen and xissimarf gen. A spectrum 
of the Fe calibration source was extracted from each CCD. The 
Mn Ka line was fit with a Gaussian profile to measure the peak 
energy. The energies were within 0.2 per cent of the expected value 
(5.895 keV) and consistent with current accuracy of the XIS energy 
scaleQ After examining for consistency, the data from the XIS-FI 
were combined to create a single spectrum. 

The PIN spectrum was extracted from the HXD data follow- 
ing standard procedures. A non-X-ray background (NXB) file cor- 
responding to the observation was obtained to generate good-time 
interval (GTI) common to the data and NXB. The data were also 
corrected for detector deadtime. The resulting PIN exposure was 
81 ks. The cosmic X-ray background (CXB) was modelled using 
the provided flat response files. The CXB and NXB background 
files were combined to create the PIN background spectrum. In 
comparison to the source plus background PIN data, Mrk 79 is de- 
tected between 12 — 25 keV. 

For the X-ray spectral analysis the source spectra are grouped 
such that each bin contains at least 20 counts (typically more). For 
display purposes data are presented in larger bins. Spectral fitting 
was performed using XSPEC vl2 . 5 . (Arnaud 1996). 

All parameters are reported in the rest frame of the source un- 
less specified otherwise. The quoted errors on the model parameters 
correspond to a 90% confidence level for one interesting parameter 
(i.e. a A% 2 = 2.7 criterion). We note the uncertainties in Section 6 
(Tabled are underestimated. Given the complexity and slow-fitting 
nature in the presented model, the uncertainties on any one model 
parameter are calculated while keeping the other variable parame- 
ters fixed. A value for the Galactic column density toward Mrk 79 
of 5.27 x 10 20 cm* 2 (Kalberla et al. 2005) is adopted in all of the 
spectral fits. Fluxes for all model components are estimated with 
the XSPEC model cf lux. K-corrected luminosities are calculated 
using a Hubble constant of Hq=70 km s _1 Mpc" 1 and a standard 
flat cosmology with M = 0.3 and S1 A = 0.7. The HXD/PIN-XIS 
cross normalisation is fixed at 1.16. 



3 CHARACTERISATION OF THE DATA SETS 
3.1 X-ray flux and spectral variability 

The seven pointed observations were conducted over ~ 7.5 years 
and span a range of ~ 3 in broadband (0.5 — lOkeV) flux. Be- 
tween 2003 and 2008, Mrk 79 underwent monitoring on hourly to 
daily intervals with RXTE (Breedt et al. 2009). The pointed obser- 
vations are compared with the RXTE light curve during this cam- 
paign (provided courtesy of P. Uttley) to determine what flux state 
is represented by each observation (FigureQJ. XMM2, XMM6, and 
SUZ1 are representative of low-flux states, with XMM6 compara- 
ble to the lowest states observed with RXTE. XMM3 is comparable 
to a high-flux state though RXTE did catch the source about 50 
per cent brighter during occasional flaring events. The pointed ob- 
servations here appear to adequately scrutinise the low-, high- and 
intermediate-flux states of Mrk 79. 

The spectra from the pointed observations are compared in 
Figure [2] The spectra appear to harden with decreasing flux and 

1 http://www.astro.isas.jaxa/suzaku/doc/suzaku_td/ 
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Table 1. Mrk 79 observation log. Each data set will be referred to in the paper by the identifier in column 1. The observatory used to obtain the data and the 
observation id are given in column 2 and 3, respectively. The start date and duration of the observation are shown in columns 4 and 5, respectively. Column 
6 shows the good time interval for each observation. The Suzaku exposure is combined for all three XIS. Column 7 provides the total source counts in the 
0.3-12 keV band (0.5 - 10 keV for Suzaku). 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Observation 


Observatory 


Observation ID 


Start Date 


Duration 


Exposure 


Counts 








(year.mm.dd) 


(s) 


(s) 




XMM1 


XMM-Newton 


0103860801 


2000.10.09 


6904 


1680 


20827 


XMM2 


XMM-Newton 


0103862101 


2001.04.26 


6935 


3590 


21454 


XMM3 


XMM-Newton 


0400070201 


2006.09.30 


24000 


14411 


238259 


XMM4 


XMM-Newton 


0400070301 


2006.11.01 


21500 


13993 


150349 


XMM5 


XMM-Newton 


0400070401 


2007.03.19 


21500 


13992 


156980 


SUZ1 


Suzaku 


702044010 


2007.04.03 


160017 


251127 


202872 












(XIS0+XIS1+XIS3) 




XMM6 


XMM-Newton 


0502091301 


2008.04.27 


90400 


58262 


136940 




Figure 1. Upper panel: The 3 — 10 keV long-term RXTE light curve is 
plotted in black. For comparison the 3 — lOkeV fluxes from each pointed 
observation are shown as red stars. Lower panel: The light curve from above 
is enlarged between days 3900 — 4605 when the pointed observations are 
more frequent. The pointed observations sample all the typical flux states 
of Mrk 79. 



a narrow emission feature at about 6.4 keV is seen at all epochs. 
Though the intermediate-flux states (XMM1, XMM4, and XMM5) 
are very comparable in broadband flux, there do appear to be some 
slight changes in the spectral shape. 

The light curves from each observation are shown in Figure[3] 
Flux variations on the short time scales (e.g. < 1000 s) probed 
with XMM-Newton are present, but are of rather small amplitude. 
The ~ 2-day Suzaku observation shows larger amplitude variations 
comparable to the fluctuations measured during the RXTE monitor- 
ing (Breedt et al. 2009). 

Spectral variability is tested for by measuring the fractional 
variability (F var ) at various energies (e.g. Pond et al. 2004). For 
observations that were at least 20 ks in duration the Fvar spectrum 
is shown in Figure[4] which confirms the low-amplitude variations 
seen in the light curves (Figure [3} and illustrates minimal energy- 
dependent variability during the XMM-Newton observations. Mod- 
est spectral variability is observed during the longer Suzaku obser- 



Figure 2. The spectra from all the pointed observations are shown af- 
ter correcting for the effective area of the detectors. The spectra become 
progressively harder with decreasing flux. The colour designation is as 
follows: XMM1 (black); XMM2 (red); XMM3 (green); XMM4 (blue); 
XMM5(violet); SUZ1 (maroon); and XMM6 (orange). 

vation as expected since the variability increases with increasing 
exposure (e.g. Markowitz & Edelson 2004). The amplitude of the 
variations appears to decrease with increasing energy. Such behav- 
ior has been seen in other AGN and attributed to various effects 
(e.g. changing spectral slope or varying the relative normalisation 
of the continuum components) (e.g Gallo et al. 2007; Ponti et al. 
2006). For our purposes, the detected spectral variability indicates 
that caution must be exercised when fitting the mean Suzaku spec- 
trum. 

3.2 Ultraviolet variability 

Simultaneous UV photometry with the OM permits examination of 
the variability at UV wavelengths. The average 2500 A luminosity 
over the six XMM-Newton observations is ~ 2.23 x 10 2S ergs -1 
with fluctuations on the ±10 per cent level (Figure[5j, significantly 
lower than seen in the X-rays. 

The power law between 2500 A and 2 keV (i.e. a ox ) is mea- 
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Figure 3. The short term light curves in 100 s bins from all the pointed 
observations. All the XMM-Newton pn (0.2 — 12 keV) light curves are pre- 
sented in the top panel. Zero marks the start time of each light curve. The 
count rates and duration of each curve are on comparable scales. XMM1 
(leftmost panel) and XMM2 (second panel from left) are short observations 
(< 10 ks). In the lower panel the Suzaku light curve for the combined XIS 
in the 0.5 — lOkeV band. Again, zero marks the start of the observation. 



Figure 5. The 2500 A luminosity (upper panel) and Aa OI (lower panel) is 
plotted for each XMM-Newton observation. The red line in the upper panel 
identifies the mean UV luminosity. Based on the photometric uncertainties 
only, the error bars on logL » are comparable to the size of the data 
points. The red line in the lower panel marks the expected a ox based on 
the UV luminosity and the shaded region corresponds to the uncertainties 
in the Just et al. relation. Negative values indicate the object is X-ray weak 
compared to the UV luminosity. 
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Figure 4. The normalised rms spectrum calculated for observations of at 
least 20 ks in duration. The XMM-Newton and Suzaku light curves used are 
in 750 s and 5760 s bins, respectively. 

sured from the simultaneous UV and X-ray observations. Accord- 
ing to Just et al. (2007), based on its UV luminosity the expected 
UV-X-ray slope for Mrk 79 is a ox (L ») = -1.25 ± 0.13. 

2500 A. 

In Figure [5] Aa ox — a ox — a ox (L °) is plotted for each 
XMM-Newton observation to illustrate th°e deviation of a ox from 
the expected value for Mrk 79 based on its 2500 A luminosity at 
the time. Most of the time Mrk 79 appears to be slightly X-ray 
strong (a ox (L ° ) > 0). XMM6 is the only observation in 

2 500 A. 

which Mrk 79 is in an X-ray weak state. According to Gallo (2006), 
marked spectral complexity above 2 keV (i.e. significant deviations 
from a power law) should be present, which is consistent with Fig- 
ure^ Though the UV luminosity is variable the plunge into an 
X-ray weak state during XMM6 is driven by X-ray variability. 



3.3 Phenomenological spectral fits 

Fitting a power law to each spectrum in the 2.2 — 4.5 and 7 — 
10 keV band and extrapolating over the 0.3 — 12 keV band (0.5 — 
25 keV for Suzaku) reveals a typical type I AGN spectrum. A soft 
excess is seen below ~ 2 keV. There is a narrow emission feature 
at ~ 6.4 keV as well as a slight excess above lOkeV (top panel 
of Figure[6]l. Qualitatively, the soft excess and the emission feature 
appear stronger with decreasing flux. The Suzaku data only extend 
down to 0.5 keV and do not show the full extent of the low-energy 
excess. 

The narrow emission feature at ~ 6.4 keV is most likely at- 
tributed to Fe Ka emission from material distant from the black 
hole. The spectra above 2.2 keV were refit with a power law and 
narrow Gaussian profile revealing the residual in the lower panel 
of Figure |6] For XMM1, XMM2, XMM3, XMM4 and XMM5 the 
simple fit appears to be acceptable. For SUZl and XMM6 positive 
residuals remain redward of 6.4 keV that could be interpreted as 
broadened Fe Ka emission from the accretion disc. For SUZl, a 
second Gaussian profile with a peak energy of 6.46 ± 0.14 keV 
(6.36 ± 0.02 keV for XMM6) and width of a = 556±^|| eV 
(142+^ eV for XMM6) can characterise the residuals. We make 
note of positive residuals extending down to ~ 4.5 keV in XMM3, 
but do not consider these as sufficient evidence for an additional 
component. 

Gallo et al. (2006) successfully fit the snap-shot observations 
of Mrk 79 (XMM1 and XMM2) with the traditional phenomologi- 
cal model of a black body and power law continuum, with a high- 
energy Gaussian profile and low-energy absorption edge. Applying 
this model to the spectrum at each epoch resulted in reasonably 
good fits {xl = 0.89 - 1.13). 

The soft excess in Mrk 79 can be parameterised as a black 
body with a varying temperature in the range 100 — 140 eV, as 
typically seen in unabsorbed AGN (e.g. Gierlinski & Done 2004; 
Crummy et al. 2006). The measured temperature is not correlated 
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Figure 6. Upper panel: The residuals (data/model) remaining from fitting 
each spectrum with a power law absorbed by Galactic absorption in the 
2.2 — 4.5 and 7 — 10 keV band, and then extrapolating over 0.3 — 25 keV. 
A low-energy excess and narrow emission feature at ~ 6.4 keV are present 
in all spectra. The Suzaku data show an excess above 12keV. Colours are 
as described in Figure [2] Lower panel: The residuals of a power law plus 
narrow Gaussian profile fit to the 2.2 — 12 keV band. Excess residuals red- 
ward of 6.4 keV are seen in SUZ1 and XMM6. 



with the source luminosity in any obvious way and does not follow 
the expected L oc T 4 relation attributed to a standard accretion disc 
of constant area (Figure |7J. The expected temperature for a stan- 
dard accretion disc around Mrk 79, assuming a black hole mass of 
5.24 x 1O 7 M (Peterson et al. 2004) and L/L E dd ~ 0.09 (fol- 
lowing Kaspi et al. 2000), peaks at ~ 15 eV. Given the high tem- 
perature predicted by the model and apparent variability, a black 
body does not provide a physical description for the soft excess in 
Mrk 79. 

A narrow Gaussian profile is a rational addition evident from 
the residuals in Figure [6] The energy is consistent with neutral 
Fe Kq emission and the flux of the line is comparable within un- 
certainties at all epochs. The average line width corresponds to a 
FWHM velocity of ~ 11000 kms -1 , about twice as high as that 
measured for H/3 from the broad line region (BLR; Peterson et al. 
1998). This suggests that the dominant component of the Fe Kq 
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Figure 7. The measured black body luminosity and temperature are plotted 
for each spectrum. The red curve is the fitted L oc T 4 relation (leaving only 
the normalisation free), which clearly does not describe the data. 



emission in Mrk 79 is not coming from a distant torus, but is more 
likely emission from the accretion disc and the inner BLR. 

The best fit absorption edge energy at most epochs is ~ 
720 eV (673 ± 6 eV for XMM6) and the depth of the feature dif- 
fers at each epoch. The energy is lower than expected from O VII 
(739 eV) and not likely associated with absorption by hot local gas 
(McKernan et al. al. 2004) since the implied velocities would be 
higher than the recessional velocity of the AGN. A more physical 
treatment is required to interpret this component of the spectrum 
(see Section[6]l. 



4 THE RGS SPECTRUM 

The RGS data were reduced using SASv9.0 (Gabrie l et al. 2003b 
and the most updated calibration files available at the date the re- 
duction was performed (October 2009). The data reduction fol- 
lowed standard procedures as in, e.g. Guainazzi & Bianchi (2007). 
The total number of background subtracted counts in the nominal 
RGS energy bandpass (0.3-2 keV) is 674 ± 19 and 727 ± 18 in the 
RGS1 and RGS2, respectively. 

Visual inspection of the RGS spectrum does not show the high 
equivalent width lines typically observed in the high-resolution 
spectra of highly obscured AGN (cf. Kinkhabwala et al. 2002, and 
references therein). The statistics of the data does not warrant de- 
tailed spectral fitting with physical models. We have therefore lim- 
ited ourselves to phenomenological simultaneous fits of the com- 
bined time-averaged RGS1 and RGS2 spectra. 

First we have performed a blind search for absorption lines. 
The spectra were fit in 100-channels intervals with a power-law 
continuum, photoelectrically absorbed by a column density held 
fixed to the value due to intervening gas in our Galaxy. At the 90 
per cent confidence level for one interesting parameter only one 
line is detected, with centroid energy E c = 698.2 ± 1.2 eV, and 
intensity / = -2.5^2' 4 x 10 -5 ctrT 1 s -2 . The difference in Cash- 
statistics between a fit with the continuum only and the fit after the 
inclusion of this absorption line is AC = 26.1. The AC corre- 
sponding to the quoted confidence level is calculated by dividing 
the distribution percentage point probability by the number of tri- 



6 L. C. Gallo et al. 



Table 2. Intensities (in units of 10 -5 cm -2 s — 1 ) and values of the R, L, 
and G diagnostic parameters for the N VI and the O VII triplets in the time- 
averaged Mrk 79 RGS spectra 



Species 1(f) I(i) I(r) R L G 

NVI 2.5 ±1.5 < 2.5 1.4 ±1.3 > 1.0 < 3.8 < 4.6 
OVII 4.9 ±1.4 < 1.8 < 0.7 > 2.1 



als, in this case the number of RGS spectral bins (3400), under the 
approximation that they are independent (this is largely valid for 
high-resolution cameras) (Svoboda et al. 2010). There is no obvi- 
ous identification for this feature. The nominal closest transition 
according to the ATOMDB database would be a resonant transition 
between the ground state and the ls4p *Pi electronic configura- 
tion of O VII. Given the lack of other detected transitions which 
could provide information on the velocity shift of the absorbing 
gas, as well as the modest statistical significance of this detections, 
we consider it as tentative only. We therefore do not discuss it fur- 
ther. 

On the other hand, several emission lines are detected once 
a similar search procedure as described in Guainazzi & Bianchi 
(2007) is applied to the combined time-averaged RGS spectra of 
Mrk 79. Most of of them are associated to the N VI and O VII 
triplets. The centroid energies of the triplet forbidden components 
are 431 ±3 and 5 74 ±4 eV, respectively. The intensity of the triplet 
components are reported in Table [2] together with the values of 
the R, L, and G diagnostic parameters (Porquet & Dub au 2000| 
IPorter & Ferland 2007t . The OVII is also detected: E c = 
67lt 2 eV, I = 1.7 ± 0.8 x 10" 5 cm -1 s -2 , while the detec- 
tion of the OviII Ly-a is only marginal (i.e. lower than the 90% 
confidence level). 

The values of the R and G parameters, as well as the high in- 
tensity ratio between the forbidden component of the O VII triplet 
and the O VIII Ly-a (4.4 ± 3.3; Guainazzi et al. 2009) indi- 
cates photoionisation as the most likely source of ionisation for 
the gas responsible for the emission features. The properties of 
the photoionised plasma are, however, only poorly constrained: 
n e < 10 10 cm" 3 , T < 10 5 K, N H > 10 20 cm" 2 . These val- 
ues are consistent with those typically measured from large-scale 
photoionised plasmas in obscured AGN (Kinkhabwala et al. 2002; 
Sako et al. 2000; Armentrout et al. 2007; Awaki et al. 2008). 



5 THE 3-12 keV BAND 

5.1 The existence of a broad line 

Figure [6] shows excess residuals in XMM6 and SUZ1 redward 
of 6.4 keV that could be described by relativistically broadened 
iron emission. When fitting broadband continuum models such as 
ref lionx (Ross & Fabian 2005) to the 0.3 - 12 keV range, the 
fit will be driven by the soft excess since the statistics are high- 
est below 2 keV. Thus the reflection and blurring model param- 
eters depend on the interpretation for the soft excess. Here, only 
the spectrum above 3 keV is examined for the existence of broad- 
ened iron emission. Furthermore, only SUZ1 and XMM6 are con- 
sidered as residuals that can be described with a broad Gaussian 
profile (Section [3~3l > are clearly seen at those epochs. This enhances 
the potential to distinguish the power law continuum and reflection 
component. 

Initially the power law component is replaced with pexrav 



Table 3. Parameter values for relativistic broad line (Laor profile) fits to 
SUZ1 and XMM6. Values that are linked between epochs appear in only 
one column. Flux units are erg cm -2 s -1 . The / superscript identifies 
fixed parameters. 



Parameter SUZ1 XMM6 



E rest (keV) 


6.77 ±0.12 


6.47 ±0.10 


a 


3' 




Rin ( ^ g) 


6.23i^ 




Rout (r g ) 


400 / 




in 


21+_i 




l°g-Fb.5-10 keV 


-12.63tnf> 


19 ^9+0.05 



(Magdziarz & Zdziarski 1995) to describe the reflection of an ex- 
ponential cut-off power law from a neutral disc. The model does 
not include an emission line, which is inserted separately. Accord- 
ing to the best-fit equivalent width (EW) measurements of the 
narrow Fe Kq line at each epoch, the reflection fraction (R = 
EW/180eV; George & Fabian 1991) is set to 0.58 and 0.89 for 
SUZ1 and XMM6, respectively. In addition, two narrow Gaussian 
profiles representing Fe Kq and Fe Kf3 are included at fixed ener- 
gies and relative normalisations (Fe K/3 = 0. 17Fe Ka). The model 
is statistically acceptable (x 2 /dof = 0.950/2007) and fits the nar- 
row core seen in Figure [6] (upper panel), but leaves broad residuals 
at both epochs (top panel Figure [8). Two additional Gaussian pro- 
files corresponding to Fe XXV and Fe XXVI are included, but do not 
improve the fit or residuals significantly (middle panel Figure[8). In 
addition, the flux of Fe XXVI is consistent with zero in XMM6. 

Replacing the two high-ionisation lines with a single Laor pro- 
file (Laor 1991) improves the residuals significantly (lower panel 
Figure [8} . The number of free parameters are the same as in the 
model with high-ionisation lines since only the energy and normal- 
isation are free to vary. The best-fit energies are 6.60 ± 0.08 keV 
(SUZ1) and 6.35 ± 0.04 keV (XMM6). 

In a second fitting attempt with the Laor model, the disc in- 
clination (i) and inner disc radius (Ri n ) are allowed to vary, and 
the emissivity index (a) remains fixed. This was done in order to 
estimate some of the parameters for when more complicated mod- 
els are attempted (see Section |6). The derived fit parameters are 
shown in Table [3] The differing line energies at each epoch could 
be suggestive of a change in disc ionisation. These preliminary 
fits argue against a maximally rotating black hole in Mrk 79 since 
Rin > 1.2 r g (where lr s = GM/c 2 ). 

Alternatively the 3 — 12 keV spectra are fitted using a power 
law for the intrinsic continuum and ref lionx for the reflection 
component. Blurring of the spectrum, which is thought to originate 
close to the black hole, is accomplished with kerrconv (Bren- 
neman & Reynolds 2006). In this case, the blurring and reflection 
parameters: inclination (i), iron abundance (Ape), black hole spin 
(a), outer disc radius (Rout) are not expected to vary significantly 
between the two observations hence are linked. The power law con- 
tinuum (r and normalisation) as well as the ionisation parameter 
and normalisation of the reflector are permitted to vary. The qual- 
ity of fit is equivalent to the Laor model (\l /dof = 0.836/2005), 
but provides a more realistic interpretation for the data. The mea- 
sured spin parameter is a — 0.7 ± 0.1. In Figure|9]variations in the 
quality of 3 — 12 keV band fit associated with the spin parameter 
are plotted for the model described. Moderately high spin values 
(a ~ 0.40 — 0.85) are favoured. We also examine the influence of 
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Figure 8. The residuals (data/model) remaining from fitting the XMM6 
(orange) and SUZ1 (maroon) spectra between 3 — 12keV with a pexrav 
continuum and line models as indicated in each panel (see text for details). 



Figure 9. The quality of the fit is tested against the black hole spin pa- 
rameter showing that extreme (high and low) spin values are discriminated 
against. The minimum Ax 2 is at a ~ 0.7. The red dashed line marks the 
3<r level. 
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Figure 10. The Ax 2 = 2.30, 4.61, 9.21 contours between iron abundance 
and black hole spin. 



iron abundance on the spin parameter. In Figure [10] the uncertain- 
ties in iron abundance and spin value are plotted together showing 
that they influence each other rather modestly. 



5.2 Consideration of absorption models 

Miller et al. (2007, 2009) have described AGN X-ray spectra as be- 
ing shaped by a composition of ionised absorbers partially covering 
the primary source. In doing so, their model diminishes the neces- 
sity of a very broad relativistic iron lines in some cases. The addi- 
tion of either one or two ionised absorbers to modify the continuum 
of Mrk 79 produces an equally good fit as the models described 
above. The appearance of the emission feature is reduced, but not 
dismissed (Figure [TTJ. The addition of a Gaussian or Laor profile 
improves the residuals. In terms of a Gaussian profile the best-fit 
parameters for the SUZ1 and XMM6 data are respectively: E — 

6.54 ± 0.08 and 6.34 ± 0.04 keV; a = 171±H and 1991^ eV; 

l 



and flux F = 1.23 ±0.04 and 2.93 ±0.03 x 10" 13 erg cm" 2 s 
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Figure 11. The residuals (data/model) remaining from fitting the XMM6 
(orange) and SUZ1 (maroon) spectra between 3 — 12 keV with ionised par- 
tial covering continuum and line models as indicated (see text for details). 

6 A MULTI-EPOCH SELF-CONSISTENT MODEL FOR 
Mrk79 

The above analysis demonstrates that the spectra of Mrk 79 are 
rather complex and a composition of multiple effects. In this sec- 
tion we will treat each phenomenon with more physical models 
while attempting to fit the data from all epochs together. The broad- 
band, multi-epoch analysis is driven by the findings in the previous 
sections. 



6.1 The baseline model: ionised disc reflection 

The positive residuals above a power law continuum between 6 — 
7keV seen in the low-flux states can be attributed to broadened 
Fe Ka emission from the accretion disc near the black hole. This is 
typically the strongest signature of reflection in AGN X-ray spectra 
(e.g. Miller 2007). Partial covering models can also describe the 
continuum in the low-flux states, but do not eliminate the residuals, 
and consequently require a second component. 

Illumination of a modestly ionised accretion disc by a power 
law continuum source will also generate fluorescence of low- 
ionisation species (e.g. Fe La; Fabian et al. 2009; Zoghbi et al. 
2010) below ~ 2keV (Ross & Fabian 2005; Ballantyne et al. 
2001). Once blurred for relativistic effects close to the black hole 
it has been shown that the low-energy reflection component can 
mimic the black body-like soft excess seen in most type I AGN 
(e.g. Crummy et al. 2006). 

As a starting point for fitting all the Mrk 79 spectra, the power 
law plus ref lionx model used to fit the 3 — 12keV band dur- 
ing the low-flux states (i.e. XMM6 and SUZ1, see Section [5} is 
applied to the broadband spectrum at all flux states. Even though a 
prominent broad line is not obvious in the high-flux state spectra, 
expectations are that a reflection component is still present (i.e. the 
power law source still illuminates the accretion disc), but that its 
flux is diminished relative to the continuum. As with the analysis 
in Section|5] kerrconv is used to blur the reflection component. 
One advantage of fitting multi-epoch data simultaneously is that 
parameters not expected to vary on the time scales probed can be 
linked. In this case, the reflector and blurring parameters: inclina- 



tion (i), iron abundance (Ap e ), black hole spin (a), outer disc radius 
(Rout) are not expected to vary on time scales of 8-years hence are 
linked between all epochs. 

At all epochs, the inner disc radius (Ri„) is assumed to be lo- 
cated at the innermost stable circular orbit (ISCO), hence fixed at 
Rin =1.0 nsco- Initial fits always place the outer disc radius at 
large distances so for simplicity R ou t is fixed at 100.0 nsco- The 
emissivity profile of the disc is a power law in radius (J(r) oc r~ a ) 
where a = 3 is commensurate with a simple lamp post model. The 
index is left free to vary at each epoch. Including a different in- 
dex for the inner and outer portion of the disc (i.e. a broken power 
law emissivity profile) was initially considered, but deemed unnec- 
essary. The two indices were usually comparable and/or the break 
radius separating the inner and outer disc regions was very large 
(e.g.> 150 r g ). The spin parameter is fixed to a — 0.7, the value 
calculated when fitting only the low-flux states (Section [5}- Like- 
wise the iron abundance is set to 1.2 x. This was primarily done 
to keep from generating arbitrarily high spin values that would 
artificial smooth the spectrum at low-energies (i.e. the soft ex- 
cess) where statistics are best. The ionisation parameter of the disc 
(f; = L x /nr 2 , where n is the hydrogen number density), power law 
continuum (F and normalisation), and normalisation of the reflec- 
tor are permitted to vary at each epoch. 

The blurred reflection and power law serves as a reasonable 
baseline model for all the spectra (xS/dof = 1.247/7425). It re- 
sults in similar residuals as found with a black body and power law, 
specifically a narrow emission feature at Fe Ka energies and sharp 
wiggles at lower energies. The residuals in the Fe Ka complex are 
treated by adding two narrow Gaussian profiles to mimic emission. 
The energy of one Gaussian profile is fixed at 7.058 keV represent- 
ing emission from Fe K/3. The second Gaussian profile, likely due 
to distant Fe Ka emission, is free to vary in energy, width and nor- 
malisation. The width of the Fe K/3 profile is linked to the free line. 
In addition, the flux of the Fe K/3 line is constrained to be 0.17 
times that of Fe Ka. Narrow emission lines from highly ionised 
species of iron (Fe XXV and Fe XXVI) were considered unneces- 
sary in Section [5] and consequently not modelled here. However 
the residuals of each spectrum are examined for possible features. 
The normalisation of the Fe Ka line is linked between all the spec- 
tra as there was no indication of it being variable from one epoch 
to another (see Section |33t . In total, the addition of two Gaussian 
profiles account for an extra three parameters in the fit. 

6.2 The photoionised emitter and warm absorber 

The RGS analysis in Section [4] reveals possible absorption, but 
predominately emission, from an ionised plasma. The domi- 
nance of the forbidden component in the O VII triplet and high 
O VIl(f)/0 VIII Lya point to photoionisation as the primary source 
of ionisation. The XSTAR (Kallman et al. 1996) analytical models 
warmabs and photemis are used to model the warm absorption 
and thermal (i.e. recombination and collisional excitation) emis- 
sion, respectively. The ionised plasma is assumed to have a density 
of 10 4 cm -3 consistent with the density estimated in the RGS anal- 
ysis. Element abundances were frozen at solar values. The absorber 
is assumed to cover the power law and reflection components. 

There is significant freedom in how the ionised plasma is mod- 
elled at each epoch and different approaches are considered here. 
For starters, the parameters of the photoionised emitter are linked 
between all epochs. One possibility to handle the absorber is to 
link the column density at all epochs and allow only the ionisa- 
tion to vary. This is a considerable improvement over the base- 
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line model discussed above (xl/dof = 1.028/7412). Statistically, 
a better quality of fit is achieved when the ionisation of the ab- 
sorber is linked and the column density is free to vary at each 
epoch (Ax 2 = ^84 for the same number of variables). How- 
ever, in both cases sharp residuals remain at low energies in some 
spectra. Allowing both the column density and ionisation to vary at 
each epoch is a significant improvement over fits where only one 
of the parameters is free to vary (Ax 2 = 78 for 7 additional free 
parameters). In addition the residuals at low energies are also im- 
proved (x 2 /dof = 1.006/7406). Allowing the photoionised emit- 
ter to vary in ionisation and normalisation at each epoch is only a 
marginal improvement to the fit (x„/dof — 0.999/7394). How- 
ever, it is physically feasible for the emitter to vary on such time 
scales and since the absorber is permitted to vary this is a consis- 
tent approach. 

The described model produces an acceptable fit to the multi- 
epoch spectra of Mrk 79 (xl/dof = 0.999/7394). In Figure [HI 
the model used to fit the intermediate-flux state of Mrk 79 (XMM1) 
is shown as well as the residuals resulting from the fit for each 
spectrum. XMM3, XMM4, and XMM5 benefit from simultaneous 
RXTE observations in which a 10 — 20keV flux could be mea- 
sured. These points are included on Figure [T2] (blue diamonds) to 
show the very good agreement between the model predicted flux 
and the measured flux, but they were not fitted. The model param- 
eters at each epoch are presented in Table [4] We reiterate that pa- 
rameter uncertainties will be underestimated in this section. Given 
the complexity of the model, when calculating uncertainties for a 
given parameter the other variables were held fixed. Consequently, 
the variability in a given parameter from one epoch to the next is 
not as significant as it appears, and examining general trends in the 
variations is perhaps more valuable. 

A complexity that arises in accepting the parameter uncertain- 
ties at face value is that the model predicts an apparent correlation 
between power law flux and photon index, in the sense that the in- 
trinsic spectra are harder in the low-flux states. This prediction is 
different from the general interpretation of reflection or absorption 
models where the intrinsic slope is constant, but the apparent slope 
changes due to varying contribution of the separate components. 
We reconsider the error estimation for T in the XMM6 and SUZ1 
spectra individually. While maintaining all multi-epoch determined 
parameters fixed (e.g. inclination, abundances) and varying the re- 
maining parameters, we estimate the more likely uncertainty in F 
at each low-flux state. In both cases, Ar ~ 0.1 at a 3<r level. While 
this alleviates the concern for significant intrinsic power law shape 
changes during SUZ1, it is still problematic during XMM6. Intrin- 
sic spectral changes in the power law component are predicted if 
the primary emission originates from a jet or Comptonisation, but 
the large changes from XMM6 to other epochs are difficult to rec- 
oncile. 



7 DISCUSSION 

7.1 The photoionised plasma 

The state of the photoionised emitter and warm absorber shows no 
correlation with the flux of the power law emitter (Figure [TJt. This 
alone does not dismiss a connection between the power law and 
plasma since the sparse data makes it impossible to track delays 
caused by light travel time effects and recombination time scales. 
Of course, the light travel time delays would be zero for the line- 
of-sight absorber. 




i 10 
Observed energies (keV) 




Observed energy (keV) 

Figure 12. Upper panel: The model used to fit the spectrum of XMM1. 
XMM1 represents and intermediate flux state for Mrk 79. The model is 
similar for all the spectra so only one spectrum is shown for clarity. Lower 
panel: Spectral residuals (in sigma) at each observation resulting from the 
simultaneous multi-epoch fit described in the text and Table [4] The blue 
diamonds in the residuals are based on simultaneous 10 — 20 keV flux 
measured with RXTE. They are shown to compare the model predicted flux 
above 10 keV with data, but were not part of the fit. 



On average the emitter exhibits a higher ionisation parame- 
ter than the absorber indicating that they are different structures. If 
both are ionised by the AGN and have similar densities, then the 
absorber is the more distant of the two. Using the ionisation param- 
eter to estimate the distance (r = (Z^/n^) 1 / 2 ) between the AGN 
and absorber places the absorber at ~ 2 x 10 19 cm (6.5 pc). The 
emitter would be approximately half as far from the AGN. Both are 
too distant from the AGN to detect correlations with AGN lumi- 
nosity from our data. 

The ionisation parameter of the emitter is comparable at all 
epochs except for the lowest flux state (XMM6) where the level of 
ionisation decreases. It is not possible to determine if this change 
is driven by some prolonged decrease in the continuum flux of the 
AGN at some point in the past. The emitter shows no significant 
changes in flux level. 

The column density of the near-neutral warm absorber is on 
the order of 10 21 cm -2 , comparable to hydrogen column densities 
seen in our own Galaxy and could be attributed to a similar struc- 
ture in Mrk 79. The column density does exhibit some variability 
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Table 4. Best-fit model parameters for the multi-epoch modelling of Mrk 79. The model component and model parameter is listed in Columns 1 and 2, 
respectively. Each subsequent column refers to a specific epoch. The inner (Rin) and outer (R, ut) disc radius are given in units of the radius of marginal 
stability (r ms ), where we have further assumed that r ms = nsco- The log of the ionisation parameter is shown for both ionised plasmas and the warm 
absorber column density is in units of log(N/10 22 cm" 2 ). Values that are linked between epochs appear in only one column. The / superscript identifies fixed 
parameters. Fluxes are corrected for Galactic absorption and are in units of erg cm -2 s —1 . 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Model Component 


Model Parameter 


XMM1 


XMM2 


XMM3 


XMM4 


XMM5 


SUZ1 


XMM6 


Power law 


r 


1.76 ± 0.01 


1.95 ± 0.01 


2.00 ± 0.01 


1.81 ± 0.01 


1.86 ± 0.01 


1.61 ± 0.01 


1.49 ± 0.01 




loJT Pn r inL-lf 
±VJ & * U.O — lUfeev 


-10.835 ± 0.007 


—10.498 ± 0.006 


-10.440 ± 0.002 


-10.526 ± 0.002 


-10.593 ± 0.003 


— 10.712 ± 0.002 


-11.164 ± 0.003 


Warm Absorber 


log rtH 


-0.68 ±0.03 


-0.83 ± 0.03 


-1.11 ± 0.02 


-1.10 ± 0.02 


-1.10 ± 0.02 


-0.86 ± 0.02 


-1.08 ±0.06 






1 03+°" 


- 1 ■ ^ - . 08 


0.98 ± 0.04 


81 + °™ 

w - ui — 0.04 


1.00 ± 0.04 


w - t " J — 0.04 


1.73 ± 0.03 


Photoemitter 


i°g? 

l°g -Fn.5-10fceV 


2.09±J;jg 
-13.03l°; 2 7 ? 


-13.001S 


2 07+ - 25 
z - u <_0.09 

-12 97+ 21 
lz ' a '-o.4l 


9 ,0+0.41 
z - ±o -0.21 

ic>.io_ 64 


-13.0912:50 


2.08+g;" 
-13.431?;^ 


, .,g+0.05 
± - OO -0.03 

-12.80 ±0.03 


Narrow Lines 


E FeKa (keV) 
E FeKfl (keV) 
a (eV) 
log F FeKa 


6.40 ± 0.01 
7.058' 

-12.54 ±0.03 














Blurring 


<\ 

a 

Rin 
Rout 

s(°) 


4 62 +0.31 
-0.27 

0.7' 
1.0' 
100' 
24 ± 1 


4 16 +0 - 51 


4.83 ±0.11 


3.78 ±0.13 


5.53 ±0.13 


3.2712:1? 


4.13 ±0.10 


Reflection 


£ (erg cms -1 ) 
A Fe (Fe/solar) 


179 ±5 
1.2' 


57 ±4 


131 ±3 


130 ± 2 


145 ±2 


177 ±6 


60 ± 1 




log F .5_ 104e v 


-11.126 ±0.014 


-11.174 ±0.029 


-10.804 ±0.005 


-11.154 ±0.008 


-10.892 ±0.005 


-11.551 ±0.017 


-11.439 ±0.006 


Total flux 


log Jb.5-l0fceV 


-10.733 ±0.005 


-10.445 ±0.005 


-10.354 ±0.001 


-10.470 ±0.002 


-10.492 ±0.002 


-10.674 ±0.002 


-10.962 ±0.002 




log F2-10keV 


-10.929 ±0.005 


-10.701 ± 0.005 


-10.631 ±0.001 


-10.687 ±0.002 


-10.736 ±0.002 


-10.834 ±0.002 


-11.083 ±0.002 



(~ 20 per cent). Such variations are not very common in type I 
AGN, but have been reported previously (e.g. Gallo et al. 2007). 



7.2 Bulk motion of the primary source and the role of light 
bending in Mrk 79 

Objects like IRAS 13224-3809 and 1H 0707-495, where light 
bending is attributed with a substantial role, often exhibit high 
values for the emissivity index (e.g. a > 6) (e.g. Ponti et al. 
2010; Zoghbi et al. 2010), indicating X-ray emission originating 
from a centrally compact region. In Mrk 79 the best-fit model was 
achieved when disc emission over a large area was allowed. This is 
evident by the relatively shallow emissivity profile measured for the 
disc. The emissivity index in Mrk 79 is typically between a ~ 3 — 5 
and this may even be slightly overestimated since, in general, limb- 
brightening effects are not considered in our blurring model (Svo- 
boda et al. 2009). The emissivity in Mrk 79 is commensurate with 
a simple lamp post model (a — 3). 

In the light bending model the primary component is normally 
perceived as an isotropic and relatively constant emitter. The ob- 
served rapid variability, as described by Miniutti & Fabian (2004), 
is driven by changes in the distance of the primary emitter from the 
black hole/accretion disc. The model makes specific predictions on 
the temporal behaviour of the reflection and power law component. 
For example, variations in the two components are correlated when 
the illuminating source is within a few r g of the black hole. At 
intermediate distances (~ 10 r g ) the reflection component remains 
rather constant while the power law can vary significantly. On long- 
time scales, as observed in this work, the power law and reflection 
components will be largely correlated assuming that the illuminat- 
ing component is emitting isotropically. In addition, light bending 
does not exclude the possibility of long term intrinsic changes in 
the flux of the primary emitter. Given that these data span nearly 



eight years, long-term intrinsic changes could take place as well. In 
Figure fj4] average flux values for each observation are plotted for 
the reflection and power law components. The figure points to vari- 
ations of similar amplitude in both components. There is a common 
trend between the fluxes of the power law and reflection, but not a 
one-to-one correlation. 

More importantly, the relative contribution of the reflection 
component in Mrk 79 is small. We estimated the reflection fraction 
(R) by comparing pexrav with simulations of the spectral models 
above 10 keV at each epoch. The power law slope, abundances and 
inclination are fixed to the values derived with the fits (TableQ. The 
cutoff energy is fixed at 350 keV, the value used byreflionx.At 
all epochs the reflection fraction was between R ~ 0.24 — 0.73, 
with the highest R value occurring for XMM6. Notably, R is vari- 
able from epoch-to-epoch and is consistently less than one. While 
reflection features are present in Mrk 79, light bending does not 
appear to be a significant factor. 

The variations in R also advocate a non-standard geometry 
for the illuminating source. If the source were simply varying its 
luminosity to produce the observed variability, there should exist a 
good correlation between the reflection and power law flux along 
with a constant reflection fraction. This is not the case. Therefore 
the standard picture in which the primary source, wherever it is lo- 
cated, only varying in luminosity with its other properties remain- 
ing unchanged is challenged. 

One can imagine a situation in which a compact primary 
source is located at a few r g , hence subject to general relativis- 
tic effects, but is moving away at mildly relativistic speeds (e.g. 
Reynolds & Fabian 1997; Beloborodov 1999). Beloborodov pro- 
posed such a model to explain the low reflection fraction in the 
black hole binary Cyg X-l. Beloborodov described the nature of 
the primary source as an active corona atop the accretion disc. The 
plasma in the active regions is heated by magnetic flares, and it can 



Multi-epoch X-ray observations ofMrk 79 11 



' 1 1 1 

( 


I < 


I ' 1 

, J 




I 
) 

[ 


' 1 ' I 


m 
CD 






9 A m ( 

HI 

s s g s s 

x x K xx 

- . , i , . , i . , . i 


■ XMM2 w 
XMM3 a 

I.I.I 




Power law flux (10 erg cm s ) 



1 2 
Power law flux (10 



erg cm " s ) 



ft. £ 
bo 1 
° *t 



Z 

60 — 



-i 1 1 1 1 1 1 1 1 1 1- 



s s 
1 1 



Power law flux (10 erg cm s ) 

Figure 13. Upper panel: The logarithm of the ionisation parameter of the 
photoionised emitter (black circles) and absorber (red squares) is plotted 
against the flux of the power law component. The average ionisation pa- 
rameters for each component are also shown (straight lines). Lower panel: 
The 0.5 — 10 keV flux of the emitter (black circles) and column density of 
the absorber (red squares) is shown against power law flux. Uncertainties 
are plotted on all quantities, but may be smaller than the data point itself. 



be ejected at mildly relativistic velocities if the blob is dominated 
by e ± pairs. As a consequence of beaming the reflection fraction is 
reduced even in the presence of light bending effects. From equa- 
tion (3) of Beloborodov (1999), R = (1 - 0/2) (1 - f3cosif/{l + 
f3) 2 , the bulk velocity of the illuminating blob (/? = v/c) is related 
to R, where i is the inclination of the system. We determine /3 at 
each epoch from the estimated reflection fraction. The maximum 
velocity will occur at the lowest R (~ 0.24) and in this case is 
/3 ~ 0.32. As the source is radio-quiet one would not expect these 
ejections to escape. For the estimated /3, the escape radius from the 
5.24 x 10 7 M black hole in Mrk 79 is w 20 r g . For higher values 
of R, one can allow for smaller ft, which would result in sources 
within an even larger distance. Light bending could have a dimin- 
ished role in Mrk 79 if the source is dominated by beaming. The 
fact that Mrk 79 tends to be X-ray strong for its UV luminosity 
(FigurefSJ is consistent with the X-rays being beamed. 

Asa consistency check we used equation ( 1 0) of Beloborodov, 



Figure 14. The 0.5 — 10 keV unabsorbed flux of the reflection and power 
law components are plotted. Variability of similar amplitude is present in 
both. 



T « 2((1 - f3 2 y 1/2 (l + /3)) -0 3 , to compare the predicted pho- 
ton index with our measured values (Table|4]l. There is very good 
agreement between the measured and predicted T (within a few per 
cent) for the high-flux observations (XMM1-XMM5). For the low- 
flux observations (SUZl and XMM6) the difference is 10 — 20 per 
cent, and alone, the F(/3) relation for this active corona model can- 
not explain the low-state spectral hardening. As discussed in Sec- 
tion 6, this is more problematic for XMM6. We note that SUZl and 
XMM6 were the longest observation, and did exhibit some spectral 
variability within the observation (Figure[4j that would affect fitting 
of the average spectrum. 

From long-term optical and X-ray monitoring of Mrk 79, 
Breedt et al. (2009) find that the optical (l/-band) and X-ray are 
highly correlated with a delay of Ol^j days. Breedt et al. find that 
the optical variations could still be produced by thermal repro- 
cessing of variable X-rays if the height of the illuminating X-ray 
source is on the order of 21.5 r g . At ~ 20 r g the light travel- 
time delay between the illuminating X-ray source and the disc is 
t = 20 r g /c = 5000 s. Consequently, a delay of r ~ 5 — 10 ks 
between the power law and reflection component is possibly de- 
tectable. Despite all the high-quality data in this study, the XMM1- 
XMM5 observations last <C 20 ks and the Suzaku data include gaps 
in the light curve every ~ 5.6 ks due to its low-Earth orbit ren- 
dering them unsuitable to search for a > 5 ks lag. Only XMM6 is 
sufficiently long enough to search for such a lag, and interestingly 
the cross correlation between the 0.2 — 0.8 keV (reflection compo- 
nent) and 1.5 — 2 keV (power law component) band does suggest 
a shift in the peak to negative times (i.e. high-energy leads low- 
energy changes). However, the XMM6 light curves suffers from 
the fact that it drops smoothly with time and lacks distinct features 
to build a correlation on (Figure[3](. Any detection of a lag with this 
light curve is judged questionable. 

Coincidently, Breedt et al. (2009) find a best-fit disc inclina- 
tion of 24° from modelling the long-term multiwavelength light 
curve. We determine a similar inclination i = 24 ± 1° in our inde- 
pendent X-ray spectral analysis of the source. 

On average, the analysis of Mrk 79 is consistent with the il- 
luminating (power law) source being beamed and located at large 
distances (.> 20 r g ) from the disc. 
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7.3 Limitations in measuring the black hole spin 

Data quality and spectral models have reached such sophistication 
that the black hole spin parameter can now be measured from the 
broadened Fe Ka line in Galactic black holes (e.g. Miller et al. 
2009) and AGNs (e.g. Brenneman & Reynolds 2006; Miniutti et 
al. 2009; Schmoll et al. 2009). 

For Mrk 79, a moderately high spin value is measured (a = 
0.7 ± 0.1) by modelling only the high-energy spectrum. Under the 
assumption that the soft excess in AGNs is the reflection compo- 
nent it becomes tempting to use the high signal-to-noise data below 
3 keV to constrain black hole spin. In principle this is fine, but one 
must acknowledge how complicated this spectral region is. In ad- 
dition to uncertainties in the calibration, there may be unmodelled 
complexities in the ionised absorber/emitter and in the reflection 
spectrum itself. These complexities could strongly affect conclu- 
sions about spin. 

Figuresl9landll0lalso seem to exclude the possibility of a max- 
imally spinning black hole. In general, it is also difficult to con- 
clusively exclude extremely high spin values since our models are 
built under the assumption that the disc is truncated at the ISCO. 
One cannot exclude the possibility that a system includes a max- 
imal spinning black hole, but that the disc is truncated at a larger 
radius. The broad line in Mrk 79 is particularly weak (Figure[8j and 
could be consistent with emission originating from large radii. 



8 CONCLUSIONS 

Broadband X-ray spectra (0.3 — 25keV) of the Seyfert 1.2 galaxy 
Mrk 79 are analysed in attempt to examine long term variability. 
The data are obtained with XMM-Newton and Suzaku at seven dif- 
ferent epochs spanning nearly eight years, and cover a factor of 
about three in broadband flux between low- and high-flux states. 
The short time scale (i.e. within an observation) flux and spec- 
tral variations are typically small making Mrk 79 useful for study- 
ing the long-term spectral variability and "average" spectrum. The 
main results are as follow. 

(1) The multi-epoch, broadband continuum is consistent with a 
primary power law flux component and disk reflection, both mod- 
ified by absorption and emission from an ionized plasma. The rel- 
ative weakness of the disk reflection features in these spectra is 
consistent with a model for stellar-mass black holes wherein the 
power law continuum arises in flares that are driven vertically from 
the disk. Variations in the illuminating and reflected components do 
not suggest that light bending dominates the observed flux trends. 

(2) Under the assumption that the inner disc is truncated at nsco 
a moderately high black hole spin parameter is found (a = 0.7 ± 
0.1). 

(3) Low-density, photoionised plasma seen predominately in 
emission is observed and consistent with arising at large distances 
from the black hole. The measured parameter values are consistent 
with those typically measured from large-scale photoionised plas- 
mas in obscured AGN. 

(4) Simultaneous UV photometry was obtained for 6 of the 7 ob- 
servations. The UV are variable to a lesser degree than the X-rays. 
Measurements of a ox show that Mrk 79 tends to be in an X-ray 
bright state compared to its UV luminosity. This would be consis- 
tent with the notion that the X-rays are beamed. Mrk 79 was caught 
in an X-ray weak state during only one observation (XMM6) for 
certain. In the X-ray weak state the 2 — lOkeV spectrum appears 
harder and deviates significantly from a simple power law. 



The modest intra-observation variability exhibited by Mrk 79 
introduces the potential to obtain a high signal-to-noise average 
spectrum of the source. The multi-epoch monitoring allows for sev- 
eral model parameters to be constrained with high precision allow- 
ing for a more accurate determination of the variable components. 
Mrk 79 does not appear to be a particularly extreme active galaxy, 
but perhaps it is more representative of a typical AGN, and conse- 
quently it is important to scrutinise it. Astro-H and 1XO will have 
substantially better energy resolution and hard X-ray sensitive than 
current X-ray missions thereby improving our ability to constrain 
the warm emitter/absorber and reflection component. 
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